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If a first-order phase transition is terminated by collisions of new-phase bub- 
bles, there will exist a period of nonequilibrium between the time bubbles 



collide and the time thermal equilibrium is established. We study the behav- 



ior of the order parameter during this phase. We find that large nonthermal 
fluctuations at this stage tend to restore symmetry, i.e., the order parame- 



ter is smaller than its eventual thermal equilibrium value. We comment on 
possible consequences for electroweak baryogenesis. 
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It has long been known that symmetry may be restored at high temperature in local 
thermodynamic equilibrium (LTE) Recently it was realized that certain nonequi- 
librium (NEQ) conditions can be even more efficient for symmetry restoration 0. An 
example of such a nonequilibrium state can arise naturally after inflation in the so-called 
preheating era [§, [|]. In fact, symmetry may be restored in the NEQ state even if it 
is not restored in the LTE state formed by thermalization of the NEQ state. Detailed 
numerical studies || confirm that fluctuations of inflaton decay products is large enough 
for symmetry restoration, as well as for several other important effects, including baryo- 
genesis ||, and supersymmetry breaking JFj, and generation of a background of relic 
gravitational waves ||. 

States with properties similar to those in preheating, namely, anomalously large fluc- 
tuations and highly NEQ conditions, can arise in other situations as well. It was sug- 
gested in Ref. || that if bubble collisions produce large numbers of soft scalar particles 
carrying quantum numbers associated with a spontaneously broken symmetry, the phe- 
nomenon of (or tendency toward) symmetry restoration may occur. The basic point is 
that bubble collisions create NEQ conditions with a large number of "soft" quanta of 
average energy smaller than the equivalent LTE temperature corresponding to instanta- 
neous conversion of the bubble energy density into radiation. Since it may require several 
scattering times for the low-energy quanta to form a thermal distribution, it is rather 
reasonable to consider the NEQ period as a separate epoch. This is generally referred 
to as the 'preheating' epoch in a manner similar to the preheating phase of slow-roll 
inflation ||. 

The tendency of symmetry restoration in NEQ conditions after bubble collisions may 
be readily understood by the following (somewhat oversimplified) reasoning. Let us imag- 
ine that particles x are produced in the bubble wall collisions and are charged under some 
symmetry group, so that their mass, m x , depends upon some scalar field (the order 
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parameter of the symmetry) as = m 2 , + g</> 2 {] Here, g represents a combination of 

numerical factors and a coupling constant. As a simple example we might assume that 
the ^-dependent mass originates from a potential term of the form V x ^ = (l/2)g<f) 2 x 2 ■ 
As opposed to the large-angle scattering processes required for thermalization, forward- 
scattering processes do not alter the distribution function of the particles, but simply 
modify the dispersion relation. This is true in NEQ conditions, as well as the familiar 
LTE conditions. Forward scattering is manifest, for example, as ensemble and scalar 
background corrections to the particle masses. Since the forward-scattering rate is usu- 
ally larger than the large-angle scattering rate responsible for establishing a thermal 
distribution, the nonequilibrium ensemble and scalar background corrections are present 
before the initial distribution function relaxes to its thermal value. These considerations 
allow us to impose the dispersion relation u 2 = p 2 + m 2 (0) for NEQ conditions. 

The leading contribution of the particles created by bubble collisions to the one-loop 



effective potential of the scalar field can be shown to be AV(4>) ~ (n / ' E)m 2 J y (f)) 0, |10] 



where n and E are the number density and the energy of the x quanta, respectively. 
We may write the potential for the NEQ configuration as AV(<p) = I?neq0 2 , where 
£>neq = gn/E. In NEQ conditions, the coefficient -Bneq m &y be quite large, indeed larger 
that the corresponding equilibrium coefficient which scales like T^ H , Trh defined as the 
temperature of the universe when the thermal spectrum of radiation is first obtained. 
Therefore, the tendency of symmetry restoration may turn out to be rather independent 
of Trh- We also notice that since the energy E scales like the inverse of the bubble 
wall width A, E ~ A -1 , one can suggest that the effect of soft particles on symmetry 
restoration is stronger for thick bubble walls. 

The aim of the present paper is to investigate numerically the issue of symmetry 

Of course \ particles may coincide with the 4> particles themselves, but in this example the colliding 
bubbles are not made from the field (j>. Otherwise there can be some effect, but the original symmetry 
will not be restored. 
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restoration in bubble wall collision. We will explicitly show that at the final stage of first- 
order phase transitions when bubble collisions occur, nonthermal quanta are produced, 
and that they tend to restore symmetry. This tendency can be quantified as a shift of 
the order parameter from its equilibrium value toward smaller values. We will also 
confirm the conjecture about the dependence of the strength of symmetry restoration 
upon the bubble wall width. Finally, we will comment on the possible implications that 
our result may have for electroweak baryogenesis. 

Let us concentrate on a theory with a single scalar field <fi (the x particles of the 
above discussions must be identified with the 0) with Lagrangian density 

C = \d^dr - \m 2 <p 2 + \c<P 3 - \\<P* - Vo, (1) 

where Vo is a constant. We introduce the dimensionless variables (f = 0/0q, t = \f\(f)Qt, 
and £ = \/~X(j)QX ) where 0o will be fixed later. In the new variables the factor A0q is an 
overall multiplication factor for the Lagrangian (m = mj v^A0cb c — c/ \<fi , V = Vq/A0q) 



dtppd^ - V(<p) 



(2) 



The overall factor will not enter the equation of motion. The final step is a choice of 
a potential, which we choose such that dV/dip = </?(</? — 1) (</? — f m ). The equation of 
motion is then 

□V? + V?(V5- 1) (<p- ¥>m) = . (3) 

With this choice of dV/d<j) the extrema of the potential are transparent: it has minima 
at (p = and (p = 1 and a local maximum at tp = ip m (we thus fix the parameter (p m to 
be in the range < (p m < 1). We shall assume (p = 1 corresponds to the true vacuum, 
i.e., V(0) > V(l). Making the connection with Eq. (|I|), we conclude that = 4>q is the 
field strength in the true vacuum, and the constants entering Eq. ([!]) are m? = ip m A0q 
and c = (1 + <p m ) A0o- We shall require the absence of cosmological constant in the true 



vacuum, V(l) = 0; this gives Vb = (1 — 2<y? m )A0Q/12. Since we consider the true vacuum 
to be at ip = 1 and the false vacuum at (p — 0, we can further restrict the parameter (p m 
to be in the range < <p m < 0.5. 

Note that only one parameter, (p m , enters the equation of motion in the rescaled 
variables. This is a key point. The evolution of any initial field configuration, ip(r = 0, £), 
for fixed ip m will be the same in the rescaled variables, regardless of the coupling constant 
A. 

The initial field configuration for the problem at hand corresponds to a set of new- 
phase critical bubbles expanding in the false vacuum. Note that the evolution of a critical 
bubble is also defined by Eq. @, and consequently it is fixed when ip m is fixed. However, 
the bubble nucleation probability is a more complicated function of the other variables 
(note that nucleation became unsuppressed when (p m —>■ 0, i.e, when the potential barrier 
disappears). The nucleation probability will determine the initial separation of critical 
bubbles (in space, as well as in time). In our numerical integration we will consider 
the mean separation of bubble nucleation sites as another free parameter of the model . 
Fixing it gives one extra constraint on the set of parameters A, 0o and tp m . 

After nucleation, new phase bubbles expand and collide. After collisions, the spatial 
distribution of the magnitude of (p resembles a random superposition of many wavelength 
modes — a configuration with large field fluctuations. It is important that the system is 
classical and can be described by Eq. (|3|) from the time of bubble nucleation, through 
the time of bubble collisions and the condition of large field fluctuations. 

The random-wave configuration is quickly established after bubble collisions; essen- 
tially it is established on the time-scale of bubble collisions since there is no small param- 
eters in Eq. @. Eventually the waves interact and LTE is established. Since transform- 
ing the NEQ distribution function into an LTE distribution function involves producing 
states with small occupation number, the coupling constant A will enter the time scale 
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for the establishment of LTE. This time scale can be very long if A is small, so the NEQ 
configuration can exist for a long time. This phase has specific properties which are the 
subject of our study here. 

First, let us recall what is expected in the final LTE state. The LTE temperature 
can be found using energy conservation 



= V = A^^m^ , (4) 



which gives 



T —I 

J- LTE — — 



1/4 



5(1 - 2<p r 
7T 2 



1/4 



A 1/4 0o&, (5) 



where b is a constant of order unity and g*{T) is the number of relativistic degrees of 
freedom at temperature T. Note that Tlte approaches zero as A approaches 0. Due to 
interactions with the medium, LTE values of the model parameters, e.g., the effective 
mass, are different than vacuum values. The value of the parameters can be calculated 
as loop corrections to the action. Most important is the change of the effective mass, 
m 2 s (T) = m 2 + XT 2 / 4. At very high temperatures m 2 s (T) becomes positive, even if the 
zero-temperature value of m 2 was negative. This is a signal that broken symmetries are 
restored at high temperatures 0. 

In the model of Eq. (HD which we consider here, the symmetry can not be restored 
again after bubble collisions, but the temperature dependent contribution to the effective 
mass will be nonzero. Using Eq. we find that it scales with coupling constant as A 3//2 , 
and tends to zero as A tends to 0. Note for what follows that the temperature-dependent 
correction to the mass can be written in more general form, as m 2 s = m 2 + 3A(0 2 — (<fi) 2 )- 

Let us now find the mean value of the field p in thermal equilibrium with temperature 
given by Eq. (|5|). To leading order in the coupling constants, the equation dV e s/dcp = 
becomes 

{<p m + 3V\b 2 ) V - (1 + if m )if 2 + ^ - (1 + if m )V\b 2 /12 = , (6) 



where terms proportional to vA are reminiscent of temperature-dependent corrections 
to the effective potential rewritten in terms of our dimensionless variables. We see that 
the solution of this equation tends to (p — 1 when A —>■ 0. In other words, the mean value 
of the field ip in thermal equilibrium (established after the phase transition is completed) 
differs very little from the vacuum expectation value if the coupling constant is small. 

We can study the process of bubble collisions and subsequent chaotization by numer- 
ically integrating Eq. (||). We defined a 3-dimensional box of size I on a grid of size 128 3 
employing periodic boundary conditions. With periodic boundary conditions every bub- 
ble in the box is mirrored by its (infinitely repeating) reflections. As the bubble expands 
to fill the box, it will collide with its reflections, and there is no need to put more than 
one bubble inside the box to study bubble collisions. So we have restricted ourselves to 
an initial configuration corresponding to just one critical bubble of the true phase in the 
box. In this case the size of the box, corresponds to the mean initial separation of 
bubbles in units of l/A 1//2 0o- We integrated the equation of motion for / = 4, 8, 10 and 
12, corresponding to progressively larger bubbles at collision time. We used (p m = 0.1 
for the only parameter in the equation of motion. 

The results for the time dependence of zero mode of the field, <£> = (<p) j is presented 
in Fig. |l|, where (. . .) means the spatial average (over grid points). We see that after 
bubbles have collided (r > 16 for I = 4 and r > 40 for 1 = 8), the zero mode does 
not relax to is vacuum value, ipo = 1, but oscillates near some smaller value. We define 
(f = ((f), where bar denotes the time average over several oscillations. We find ipo ~ 0.93 
in the case I = 8 and ipo ~ 0.87 with / = 4 at r ~ 80. Note that ipo rises slightly with r, 
which is the sign of ongoing relaxation. We do not present results for I = 10 and / = 12 
since they do not differ appreciably from the case I = 8 (<£>o at Z = 12 is larger by an 
about 0.01 than the corresponding value for / = 8). 

The deviation of (p from the vacuum value is not unexpected since a random field 
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Figure 1: Time dependence of the zero-momentum mode, (po = (tp). The dotted line 
corresponds to initial bubble separation of I — 4, the solid line corresponds to / = 8. 

of classical waves is created after bubble collision, i.e., Var(yj) = (ip 2 ) — (p) 2 , is nonzero. 
The time dependence of the variance is shown in the Fig. [2]. Note again that with fixed 
initial conditions the variance does not depend upon A, i.e., it has a nonthermal origin. 

At t ~ 80, with I = 8 we have Var(y?) « 0.036 and with / = 4 we find Var(^) « 0.08. 
Again we employ time averaging over several oscillations. At small A those values are 
much larger than its equivalent LTE value Var(<^) = T£ TE /12 (see Eq. (|[)). The fact 
that Var((p) in NEQ can exceed its equivalent LTE value by many orders of magnitude 
was the main point of Ref. which studied the preheating phase after inflation, and 
of Ref. which studied conditions following bubble collisions. Our work supports the 
claim in Ref. || that NEQ phase transitions can occur in models which contain more 
degrees of freedom than the simple toy model of Eq. ([!]). 

Let us see whether we can understand the deviation of the zero mode from its vacuum 
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Figure 2: Time dependence of the variance, (tp 2 ) — (<p) 2 - The dotted line corresponds to 
initial bubble separation of I = 4, the solid line corresponds to I = 8.. 

value by the existence of a nonzero Var(^). Let us decompose the field as <p — <p + 5(p, 
and substitute this decomposition into the equation dV/dip = 0. We find in the Hartree 
approximation 



(ip m + 3(5(p 2 ))(p - (1 + (fmjPo + V^o - (! + y?m)(5y? 2 ) = 

Assuming in addition the deviation of (fa from 1 to be small, we find 

2 - 



ip = 1 



(V) 



(7) 



(8) 



1 - (p m + 3(5ip 2 } 

Using ip m = 0.1 and the values of (Sip 2 ) inferred from Fig. |2], we find (fo = 0.93 for / = 8 
and <po = 0.87 for / = 4, which are in excellent agreement with the results presented in 
Fig. |. 

We can also understand the dependence upon /, the initial bubble separation. The 
larger the initial bubble separation, the longer bubbles will expand before they collide. 
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As a bubble expands, its wall thickness decreases. Hence, colliding bubbles in the I = 8 
calculation are thinner than in the I = 4 case. Following the discussion in Ref. ||, 
we expect the average energy of the quanta created in wall collisions to scale as A -1 , 
where A is the wall thickness at collision. Since the effect of the background on the 
effective potential scales as n/E oc A, we expect the I = 4 calculation to result in a 
larger departure from the vacuum value. This expectation is confirmed by the results 
shown in the figures. 

Even though we only examine a particularly simple model, we conjecture that a de- 
viation of ifo from its thermal equilibrium value in the aftermath of bubble collisions 
may have important consequences for some applications of first-order phase transitions, 
e.g., electroweak baryogenesis. In any scenario where the baryon asymmetry is gener- 
ated during a first-order electroweak phase transition, the asymmetry is generated in the 
vicinity of bubble walls, and a strong constraint on the ratio between the vacuum expec- 
tation value of the Higgs field inside the bubble and the temperature must be imposed, 
(0(T))/T > 1 [IT]. This bound is necessary for the just created baryon asymmetry 



to survive the anomalous baryon violating interactions inside the bubble, and may be 
translated into an severe upper bound on the physical mass of the scalar Higgs parti- 
cle. Combining this bound with the LEP constraint already rules out the possibility of 
electroweak baryogenesis in the standard model of electroweak interactions, and even 
leaves little room for electroweak baryogenesis in the minimal supersymmetric exten- 
sion of the standard model |12|]. Since the rate of anomalous baryon number violating 



processes scales like exp(— (0)/T), it is clear that even a small change in the vacuum 
expectation value of the Higgs scalar field from its equilibrium value may be crucial for 
electroweak baryogenesis considerations. Our results suggest that imposing the bound 
(0(T) E q)/T > 1 may not be a sufficient condition for successful electroweak baryogene- 
sis. NEQ effects at the completion of the phase transition may reduce the expectation 
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value of the Higgs field, thus enhancing the anomalous baryon number violating rate 
with respect to its equilibrium value, making the upper bound on the Higgs mass more 
severe. Applications of our considerations to the electroweak transition may result in 
a fatal blow to many scenarios involving extensions of the standard model where the 
baryon asymmetry is generated during the electroweak phase transitions. 

The model we consider in this paper is quite simple, but it illustrates several points. 
The most important result is that NEQ conditions following bubble collisions can have 
a dramatic effect upon the effective potential. Although the model we study is too 
simple to result in symmetry restoration, the numerical results confirms the assumptions 
made in Ref. |J about the efficiency of NEQ conditions. We mentioned a possible 
direct application of our results to electroweak baryogenesis, but we believe that the 
phenomenon of NEQ effects will have other implications as well. 

The work of EWK and AR was supported in part by the Department of Energy 
and by NASA grant NAG 5-2788. IT was supported by Department of Energy grant 
DE-AC02-76ER01545 at Ohio State. 
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